Stationary Electromagnetic Fields of a Slowly Rotating 
Magnetized Neutron Star in General Relativity 

L. RezzoUa^'^, B. J. Ahmedov^'^ and J. C. Miller^'^ 



^SISSA, International School for Advanced Studies, Via Beirut 2-4, 34013 Trieste, Italy 
'^INFN, Department of Physics, University of Trieste, Via Valeria, 2 34127 Trieste, Italy 
^Institute of Nuclear Physics, Ulughbek, Tashkent 702132, Uzbekistan 
'^AS-ICTP, The Abdus Salam International Center for Theoretical Physics, Strada Costiera, 11 34014 Trieste, Italy 
^Nuclear and Astrophysics Laboratory, University of Oxford, Keble Road, Oxford 0X1 3RH 



February 7, 2008 



Abstract 

Following the general formalism presented by RezzoUa, Ahmedov and Millerl0 , we here derive ana- 
lytic solutions of the electromagnetic fields equations in the internal and external background spacetime 
of a slowly rotating highly conducting magnetized neutron star. The star is assumed to be isolated and in 
vacuum, with a dipolar magnetic field not aligned with the axis of rotation. Our results indicate that the 
electromagnetic fields of a slowly rotating neutron star are modified by general relativistic effects arising 
from both the monopolar and the dipolar parts of the gravitational field. The results presented here differ 
from the ones discussed by RezzoUa, Ahmedov and Miller'EJ mainly in that we here consider the interior 
magnetic field to be dipolar with the same radial dependence as the external one. While this assumption 
might not be a realistic one, it should be seen as the application of our formalism to a case often discussed 
in the literature. 



1 INTRODUCTION 



The rich phenomenology observed from pulsars has long motivated the study of the electromagnetic fields 
interior and exterior to a rotating neutron star. Furthermore, since neutron stars are among the most relativis- 
tic stellar objects, the study of electromagnetic fields in strongly curved spacetimes has been the subject of 
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past and recent interest. Indeed, the general relativistic modifications to the solutions of the Maxwell equa- 
tions in the curved spacedpe of a relativistic starJiave been studiedj^ince the 1960s. The initial work 
of Ginzburg and Ozernoya, Anderson and Cohen O and of PettersonM on the stationary electromagnetic 
fields in a Schwarzschild external background geometry have revealed that, in general, the spacetime cur- 
vature amplifies the magnetic fields outside compact magnetic stars. More recent work by Geppert, Page 
and ZanniasH has reconsidered this problem also in the context of non-stationary magnetic fields. 

The^eneral relativistic effects induced by rotation of the star were first investigated by Muslimov and 
TsyganO who considered the case of a slowly rotating star surrounded by vacuum and with a prescribed 
dipole magnetic field aligned with the rotation axis. Muslimov and Harding then extended this to treat 
a non-vacuum exterior More recently, RezzoUa, Ahmedov and Millertt have considered stationary and 
non-stationary solutions of the Maxwell equations in the internal and external background spacetime of 
a slowly rotating magnetized relativistic star surrounded by vacuum and with a liipole magnetic field not 
ahgned with the axis of rotation. In an independent study, Konno and KojimaEI have also considered a 
similar scenario but have limited their investigation to the stationary solution of an aligned dipole. 

In this paper we make use of the general expressions formulated by RezzoUa et al.lll to study the case, 
recently discussed in the literature 0, in which the magnetic field interior to the star is dipolar and has the 
same radial dependence as the magnetic field outside it. While this is strictly speaking incorrect since the 
interior magnetic field cannot be determined independently of the structure of the rotating star, the error 
made with this assumption is small and is compensated by the possibility of dealing with a simple but 
instructive case. 

The paper is organized as follows: in Section^ we write the general relativistic equations for the electro- 
magnetic fields in the background metric of a slowly rotating star in the form they assume in the orthonormal 
tetrads of a family of zero angular momentum observers. In Section ^ we find the stationary general rel- 
ativistic solutions to the Maxwell equations in the spacetime internal and external to the rotating star. In 
Section ^ we summarize our conclusions. 

Throughout, we use a space-like signature (—,+,+,+) and geometrical units where c = G = 1 
(However, for those expressions with an astrophysical application we make exceptions by writing the speed 
of light explicitly.). Greek indices run from to 3 and Latin indices from 1 to 3; semi-colons denote 
covariant derivatives, while commas denote partial derivatives. 



2 THE MAXWELL EQUATIONS IN A SLOWLY ROTATING SPACE- 
TIME 



Hereafter we will assume that the contribution of the electromagnetic energy density to the total mass den- 
sity is negligibly small, even for very highly magnetized stars. Under this approximation we can solve 
the general relativistic Maxwell equations on a given fixed curved background, rather than the coupled 
Einstein-Maxwell equations. In particular, we will consider the background metric as being that of a station- 
ary, axially symmetric star approximated to first order in ltie_angular velocity fl. Using a Boyer-Lindquist 



where the differential rotation Lu{r) is caused by frame-dragging and can be interpreted as the angular 
velocity of a free falling (inertial) frame. The radial dependence of u; in the region of spacetime internal to 
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the star has to be found as the solution of the differential equation 



where we have defined 
and where 



(^r''"— ^ +4— w = (2) 
dr \ dr ) dr ' 



e-(*+^' , (3) 



= n~uj , (4) 



is the angular velocity of the fluid as measured from the local free falling (inertial) frame. In the vacuum 
region of spacetime external to the star, on the other hand, w(r) is given by the simple algebraic expression 

, , _ d(f) go^ 2J 
dt r'^ 

where J — I{M, R)n is the total angular moment of metric source as measured from infinity and I{M, R) 
its momentum of inertia (see liJ for a discussion of / and its numerical calculation). Outside the star, the 
metric (|l]) is completely known and explicit expressions for the other metric functions are given by 

e^*^'-) ^ (^1 - = e-2A('-) , r > i? , (6) 

where M and R are the mass and radius of the star as measured from infinity- 
Note that the use of a vacuum Schwarzschild metric in place of (|l])Ho should not be considered 

satisfactory since the Schwarzschild metric is intrinsically inadequate for describing physical systems such 

as pulsars in which the coupling of the magnetic field and rotation plays a crucial role in the production of 

the observable electric field. 

We start our analysis by recalling the general relativistic form of the Maxwell equations 

(7) 

F"'^.^ = 47r J" , (8) 

where the four-current J" can be decomposed into a conduction current density and a proper charge 
current density PeW°' 

J"=PeU;"+r, r^«a=0, (9) 

with ly" being the components of the conductor's four-velocity. If the conduction current is carried by 
electrons^ with electrical conductivity cr. Ohm's law can then be written in its special relativistic form 

jo = (jF^pw^" , (10) 

(A general relativistic expression can be found in 

The electromagnetic field tensor Fap can be expressed thrpuglLthe electric and magnetic four-vector 
fields i?", i?" measured by an observer with four- velocity Lj' EZI 

Fap = ^UlaEp] + rjajS-iSU'B^ , (11) 



^This is a reasonable assumption if the neutron star has a temperature such that the atomic nuclei are frozen into a lattice and the 
electrons form a completely relativistic, and degenerate gas. 
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where Tr 



IS 



a/3] 



\{Tai3 ~ Tpa) and rjap-ys is the pseudo-tensorial expression for the Levi-Civita symbol 



-9 



(12) 



with g EE det\gap\ = -e2(*+^)r'^ sin^ 6 for the metric _ 

It will be useful to consider the class of "zero angular momentum observers" or ZAMOsO. These are 
observers that are locally stationary (i.e. at fixed values of r and 9) but who are "dragged" into differential 
rotation with respect to a reference frame fixed with respect to distant observers. At first order in ft they 
have four- velocity components given by 



1,0, 0,w 



(Ua) 



EE e*''') 



1,0,0,0 



(13) 



We can now rewrite the Maxwell equations in the ZAMO reference frame by projecting the electromag- 
netic vector fields onto a locally orthonormal tetrad and indicating their components with "hatted" indices: 
i.e. 

smd(r^B^)_+e^r(sm9B^) +e^rB^,^Q, (14) 



and 



(rsin6') 



(e^r sin( 



dB^ 
dB^ 



Ksinf?) B\^, 



sm U re 



sm U Lur 



'■B') ^+ioe\ (sine's' 



sin 61 (r'^E^) ^ + e^r (^sin0E^^ ^ 



(^amOB^^ 



B" 



{iur ainO) E^ ^ — (rsir 



dE^ 
~dt 



e^+^B[^ - sin 61 e*B"*j - (weVsin6l) E\ = {e^sm.j 
(e*rB^") - e^+^B\g + sin 9 {uor^E^) ^ + u:e\ (sin 9E'>^ ^ 



47re*rsin6'J^ 
„N dE^ 



(15) 
(16) 
(17) 

(18) 
(19) 



47re*+^rsin6lJ'' ,(20) 



dE't' 



e\) —- + 47re*+^rJ'^ + Aire^Lor'^ sin 61 J* 
' ot 



(21) 



Note that, apart from the general relativistic corrections produced by the static part of gravitational field 
and proportional to the metric functions A and <i>, the Maxwell equations ([l4{)-(pTl) contain new general 
relativistic terms due to the coupling between the frame-dragging effect and the electromagnetic fields. 
This coupling can then act as a source of electric and magnetic fields. 
Taking our conductor to be the star with four- velocity components 



w 



-*(r) 



1,0,0,1) 



1,0,0, 



r2 sin^ 9 



(22) 
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we can use Ohm's law ( pXj| ) to derive the following explicit components of J" in the ZAMO frame 

J^^P.^.^-^eK (23) 

r = a(E^-^-^B^] , (24) 



n f a ujr sin 9 ~\ 
J' Ie' + -^B^ , (25) 

1 1 uir sin 6 

J0 ^ aE* + -^pe ■ (26) 

To simplify the problem, hereafter we will adopt the following assumptions. Firstly, for simplicity we 
consider the case when there is no matter outside the star so that the conductivity (7 = for r > R and 
(7 7^ only in a shell with R^^^ < r < R (e.g. the neutron star crust). Finally, we require a to be uniform 
within this shell (Note that this might be incorrect in the outermost layers of the neutron star but is a rather 
good approximation in the crust as a whole.). 



3 STATIONARY SOLUTIONS 



In this Section we will look for stationary solutions of the Maxwell equations, i.e. for solutions in which 
the magnetic moment of the magnetic star does not vary in time. Note that this does not mean that the 
electromagnetic fields are independent of time. As a result of the misalignment between the magnetic 
dipole II and the angular velocity vector ft, in fact, both the magnetic and the electric fields will have a 
periodic time dependence produced by the precession of p around O. _ . 

Our strategy in searching for the solution is based on extending the Deutsch solutioncJ (i.e. the solution 
to the Maxwell equations for a misaligned rotating sphere in a Minkowski spacetime) to a curved spacetime. 
In particular, we will distinguish between an exterior vacuum solution to the Maxwell equations (for which 
fully analytic solutions can be found) and an interior non-vacuum solution. These two solutions will then 
be matched at the surface of the star. 

A considerable simplification in the problem comes from the fact that, at first order in il, the solutions 
for the electromagnetic fields do not acquire general relativistic corrections in their angular parts, and the 
corresponding Deutsch solution can be used directly. We will therefore search for a magnetic vector field 
with components 

{r,9,(j>,x,t) = F{r) [cos x cos 6* + sin x sin 6 cos X{t)] , (27) 
(r, 9, (/), X, t) = G(r) [cos x sin 6* - sin x cos 9 cos X{t)] , (28) 

B^{r,9,<P,x,t)^Hir)sinxsinX{t) . (29) 

Imposing now that = = J'^ = (the interior of the star is a perfect conductor and the exterior of the 
star is a vacuum), the radial part of the Maxwell equations (|l4|),(|l9)-(|T]) will reduce to three equations for 
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the radial eigenfunctions F{r), G{r), and H{r) 

{r^F)^^ + 2e\G = 0, (30) 

(re*77) ^ + e^+^F = , (31) 

H -G = 0. (32) 

Equations (^0|)-(^2[) show that, in the case of infinite conductivity and as far as the magnetic field is 
concerned, the use of a slow rotation metric provides no additional information with respect to a non- 
rotating static me trie 111. 

3.1 Interior Solution 



It is important to notice that the system of equations (^0[)-(|3^ combines information about the structure and 
physics of the star (through the metric functions $ and A) with information about the microphysics of the 
magnetic field (through the radial eigenfunctions F and G). As a result, and as mentioned in the Introduc- 
tion, the general relativistic solution for the interior electromagnetic fields cannot be given independently 
of a self-consistent solution of the Einstein equations for the structure of the star. In practice, to calculate 
a generic solution to (|30|)-(^^, it is necessary to start with a (realistic) equation of state and obtain a full 
solution of the relativistic star. Once the latter is known, the system of equations (|o|)-(|3^) can be solved for 
a magnetic field which is consistent with the star's structure and corresponds to a magnetic configuration 
of some astrophysical interest. In the case of a star with constant density, for example, the exact analytical 
solution of the inner magnetic field can be found if a "stiff matter" equation of state is used for the stellar 
mattertl . 

Given a dipolar exterior magnetic field, a simple model for the interior magnetic field comes from 
assuming that the latter is continuous through the surface of the star and down to some inner radius r = . 
The use of an inner radius removes the problem of suitable boundary conditions for r ^ 0, and reflects 
the basic ignorance of the properties of magnetic fields in the interior regions of neutron stars. In this case, 
the radial eigenfunctions and Gj„ are obtained through the numerical solution of a system of coupled 
ordinary differential equations Ell . An even simpler but also cruder model, adopted by Konno and KojimaB 
for an aligned dipole, assumes that the radial eigenfunctions of the interior magnetic field have the same 
functional form as the corresponding eigenfunctions in the exterior of the star 

Fjj^ = ^i-^EXT = CiF and G^^ = GiG^^^ = CiG , (33) 

where Ci is a constant to be determined through the imposition of boundary conditions. This is the model 
which we will also consider hereafter as a simple application of the generic equations derived by Rezzolla, 
Ahmedov and MillerlJ 

The form of the internal electric field is straightforward to derive in the absence of conduction currents. 
In this case, in fact. Ohm's laws (^) and ( |24| ) yield the simple expressions 

^, ^ ujTsml^g ^ e-('^+^)rsm0 _ ^ ^.^^ ^ _ sin ^ cos 9 cos A) GiG , (34) 
ce^ c 

ujf Sin Q " 6 ^ rp gj^j^ Q 
= J — = Lj (cos X cos + sin X sin 9 cos A) GiF , (35) 
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(36) 



where we have taken into account that pe = 0{ijj) and that the contribution proportional to uupe is therefore 
of higher order. Note that, apart from the red-shift correction proportional to e~*, equations (^4|)-(^6|) are 
the same as those in flat spacetime with Q, being replaced by the effective fluid velocity measured by a free 
falling observer u). 

3.2 Exterior Solution 



The exterior solution for the magnetic field is simplified by the knowledge of explicit analytic expressions 
for the metric functions $ and A. In particular, after defining = e* = e^^ = (1 — 2M / r)^/^, the system 
(^0[)-(^2[) can be written as a single, second-order ordinary differential equation for the unknown function 
F 

d 
dr 



2F = 0. 



(37) 



Let now x = 1 — r/M, then equation (BTh becomes 



_d_ 

dx 



\ + x 
I — X J dx 



[i-xyF 



2F^0, 



(38) 



which has an exact solution expressed through the Legendre functions of the second kind . In partic- 
ular, the radial eigenfunctions F{r), G{r), and H{r), are found to be 



F(r) 
G(r) = 



4M3 
37V 



IniV^ 



2M 
r 



1 



M 
r 



AA'Pr 
H{r) = G{r) . 



M 



(39) 

(40) 
(41) 



Using expressions (^9|)-(|4l|) we can now determine the value of the matching constant Ci by requiring that 
the radial magnetic field is continuous across the surface of the star, i.e. that [B''~{r = R)]jj^ — [B'^{r = 



R)]ext- ^ result, we obtain 



1 



(42) 



Collecting all of the expressions for the radial eigenfunctions, we obtain the following expressions for 
the components of the stationary magnetic field 



B" = 



B^ 



4A73 
3iV 

3A^ 
4Af2r 



IniV^ 



2M 
r 



— \nN^ + ^ + 1 



M 



M 



iV2 
1 



(cos X cos + sin X sin cos A) ^ , (43) 
(cosxsin^ — sinxcos^cos A)^ , (44) 
(sin X sin A) ^ . (45) 
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The search for the form of the electric field is much more involved than for the magnetic field. Using 
as a reference the Deutsch solution, we look for the simplest solutions of the vacuum Maxwell equations in 
the form 



= (/i + /s) [cos x(3 cos'^ 6 — 1) +3 sin x cos A sin 6 cos 0] 



(46) 



e_{h + h) 



cos X sin 20 — (l — 2 sin 6) sin x cos A 



(/2 + h) - [gi + 92) 



sin X cos 6 sin A 



(51 +52) sin X cos A , (47) 



(48) 



where the unknown eigenfunctions /i — /4, and gi — 172 can be found as solutions to the vacuum Maxwell 
equations and have radial denendence only. Skipping here the details of the derivation of the expressions 
for the radial eigenfunctions (whose explicit form is presented in the Appendix) we here report only the 
final form of the components of the stationary vacuum electric field external to the misaligned magnetized 
relativistic star 



E'' ^ 



1 



(S^^' + ^^2) (2r - 3M) In + {Qt^ - 3Mr + M^)] 



\nN^ + 3M*ujr^ \ [cosx(3 cos^ 9 — 1 ) + 3 sin x cos A sin 9 cos 9] ^ 



(49) 



E" 



C3 



[3M*r3 (r^ - iMr + 2M'^) InTV^ + 2M^r'^ (Sr^ - 6Mr + M^)] 
X [2 cos X sin 9 cos 9 — (cos^ 9 — sin^ (?) sin x cos A] ^ 



2M 



M 
1 + — 
r 



(sinxcos A)/i 



2M^r'^Nc 



( 



7C0 



[3M*r^ (r^ - 3Mr + 2M'^) \nN^ + 2M'^r^ (Sr^ - 6Mr + M^)] 
3(Dr 



2M 
r 



1 



■(sinxcosSsin A)/i 



(50) 



(51) 



where C2 and C3 are arbitrary constants to be determined through the imposition of boundary conditions. 

Expressions (|4^)-(pl]) confirm that the general relativistic dragging of reference frames introduces a 
new contribution to the form of the electric field which does not have a flat spacetime analogue. This effect 
is 0{uj) and therefore is present already in our first-order slow rotation approximation. This is in contrast 
with what happens for the magnetic fields, where higher order approximations of the form of the metric are 
necessary for frame dragging corrections to appear 
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The values of C2 and C3 can be found after imposing the continuity of the tangential electric field 
across the stellar surface. Using (^4|)-(^6|) as solutions for the internal electric field and imposing that 

[E%r = R)U = [E'{r = R)],,^, as well as [E^r = = [E^r = R)],,^, yields 

^ ^ 9i?^ [\nNl +2MlR{\ + M/R)\ 

^ 8Af 3 (3i?2 _ 67\fi? + M2) + 12M2i? (i?2 _ 3Mi? + 2Af 2) In Nl/ 

2AP [inNl +2M/R] 
^ ^30Mi?2 _ 60M2i? + 10M3 + (15i?3 - 45Mi?2 + SOM^R) InN^ ' "^^^^ 

with iV2 = 7V2(-j. ^ ^) ^ 1 _ 2M/i?. 

Note that in the limit of an aligned dipole (x = 0), equations (E)-(^o|) reduce, after some algebraic 
manipulations, to the vacuum solutions found by Konno and KojimaB . 

E^ = [6M*r3 (2r - 3M) ln7V2 + 4AfV (6r2 - 3Mr + M^)] 



2 



In N-" + SArajr-" } (3 cos^ 0~l)^i, (54) 



E^ = ^-^5 1 C \3M\^ (r2 - 3Mr + 2M2) lniV2 + 2MV2 (3r2 - 6Mr + M^)] 



sin 6 cos 6fi , (55) 
z I 

whereJjhe constant C is related to the integration constant C2 given by formula (2. 11) of Konno and Kojima's 
paperB by 

f5ujr^ Ain \ C2 ojr^ 



4 CONCLUSION 



In a recent related workB , we have presented analytic general relativistic expressions for the stationary 
electromagnetic fields in the crust and vacuum region of a slowly-rotating magnetized neutron star There, 
the star was considered isolated and in vacuum, but no assumption was made about the orientation of 
the dipolar magnetic field with respect to the rotation axis. As an application of the general formalism 
developed, we then found an exact solution for the interior magnetic field for a star with uniform density 
and a stiff matter equation of state. 

In this paper we have instead considered the case of an internal magnetic field with a radial dependence 
which is the same as that of the external dipolar magnetic field. Strictly speaking this cannot be entirely 
correct since the interior magnetic field cannot be determined independently of the structure of the rotating 
star. However, this remains an interesting application of jOur equations to a case which has a Newtonian 
analogue and has recently been discussed in the literature^ . 

The solutions presented in this paper provide a lowest order analytic form for the electromagnetic fields 
in the spacetime of a slowly rotating misaligned dipole and agree, in the case of an aligned dipole, with the 
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results of Konno and KojimaBl. While our solutions have been found with some simphfying assumptions, 
they also allow the major features of a realistic solution to be seen and could therefore be used in a variety 
of astrophysical situations. 
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APPENDIX 



For completeness, we here report the explicit forms of the functions /i — f^, and gi , 32 used in expres- 
sions (p6|)-(p8|) for the components of the external electric field 



/i 

/2 
/3 

u 



92 



n 



C2 



2r 



2Af2 2M 



-4 



C.N 



1 ]\nN^ -2 ^-^ 



2r \ , , 2AP 2M 



~ SAPcN 



2M 



InN^ + — ( 1 

r \ r 



M 



3ujr 
' SM^cN 



where + 52 = {uj/^)gi. 



r 

2AP 



\ , . 2M / M 

\nN^ H IH 

r V r 



2M2 



(57) 
(58) 

>f^, (59) 
(60) 
(61) 
(62) 
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